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AN ANALYSIS OF THE SECONDARY STRUCTURE OF THE 
MITOCHONDRIAL LARGE SUBUNIT rRNA GENE (16S) IN 
SPIDERS AND ITS IMPLICATIONS FOR 
PHYLOGENETIC RECONSTRUCTION 
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ABSTRACT. We investigated the pattern of molecular variation with respect to secondary structure in 
the 16S ribosomal RNA gene and its phylogenetic implications for arachnids with a focus on spiders. 
Based on a model by Gutell et al. (1996), secondary structures were proposed for the 3' half of 16S in 
the mygalomorph spider Aptostichus atomarius. Models were also constructed for a hypervariable length 
of the 16S in three other arachnids, which revealed a trend of stem and loop reduction in more advanced 
arachnids. Using a simple statistical approach to compare functional regions, we found that internal and 
external loops are more variable than stems or connection regions. Down-weighting or excluding regions 
which code for the more variable loops improved tree topologies by restoring the monophyly of the genus 
Aptostichus , a group supported by combined 16S, COI, and morphological data in other analyses. This 
study demonstrated the utility of considering secondary structure for DNA sequence alignment and phy¬ 
logenetic reconstruction in spiders. 
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The most important aspect of any molecular 
phylogenetic study is, unequivocally, gene 
choice. Choosing a gene that strikes the ap¬ 
propriate balance between molecular conser¬ 
vation and variability is essential to the suc¬ 
cess of phylogenetic reconstruction. Without 
question the functional role of a gene and its 
rate of evolution are tightly coupled although 
within genes, this relationship is variable be¬ 
cause regions of a single gene can evolve at 
different rates. This tight relationship between 
rate of nucleotide substitution and gene com¬ 
ponent “form and function” results in a single 
gene being useful at different, sometimes 
quite disparate phylogenetic levels. Such mul¬ 
tilevel phylogenetic utility, strongly tied to 
gene function (i.e., the functional secondary 
structure) is particularly true for the mito¬ 
chondrial 16S rRNA gene in arthropods 
(Flook & Rowell 1995), the gene of focus for 
this paper. 

1 Current address: Stacey I). Smith, Department of 
Botany, 430 Lincoln Dr., University of Wisconsin, 
Madison, WI 53706-1481. E-mail: sdsmith4@ 
wisc.edu 


The 16S rRNA gene, which encodes the 
mitochondrial large ribosomal subunit (mt 
LSU) in animals, has been employed exten¬ 
sively to explore phylogenetic relationships in 
arthropods at most phylogenetic levels [e.g., 
ordinal (e.g., Flook & Rowell 1995), familial 
level (e.g.. Black & Piesman 1994) and the 
genus level and below (e.g., DeSalle et al. 
1992; Bond et al. 2001)J. The wide range in 
utility of 16S at various taxonomic levels sug¬ 
gests that the differential rates of molecular 
evolution within 16S, due to varying func¬ 
tional constraints, greatly affect its phyloge¬ 
netic utility. Thus, understanding the function¬ 
al roles of different portions of the gene 
through secondary structure modeling should 
lead to a more robust use of 16S in phyloge¬ 
netic studies. 

The main focus of this paper is the second¬ 
ary structure of 16S in spiders and other 
arachnids. One of the first models of 16S sec¬ 
ondary structure in arthropods was created by 
Clary & Wolstenholme (1985) using the Dro¬ 
sophila yakuba sequence. Subsequent models 
by Gutell and colleagues (Gutell & Fox 1988; 


44 



SMITH & BOND-SPIDER 16S SECONDARY STRUCTURE 


45 


Table 1.—Spider taxa sampled. 


Higher level 
classification 

Species 

Origin 

Reference/GenBank 
accession # 

Mesothelae 

Liphistiidae 

Heptathela nishihirai Haupt 1979 

Japan 

Huber et al. 

1993 

Opisthothelae 

Araneomorphae 

Clubionidae 

Clubiona pallidula (Clerck 1757) 

Austria 

Huber et al. 

1993 

Ctenidae 

Cupiennius coccineus F. Pickard-Cam- 
bridge 1901 

Costa Rica 

Huber et al. 

1993 

Ctenidae 

Cupiennius getazi Simon 1891 

Costa Rica 

Huber et al. 

1993 

Ctenidae 

Cupiennius salei (Keyserling 1877) 

Mexico 

Huber et al. 

1993 

Ctenidae 

Phoneutria boliviensis (F.O. Pickard- 
Cambridge 1897) 

Costa Rica 

Huber et al. 

1993 

Lycosidae 

Pardosa agrestis (Westring 1861) 

Austria 

Huber et al. 

1993 

Pisauridae 

Dolomedes fimbriatus (Clerck 1757) 

Austria 

Huber et al. 

1993 

Pisauridae 

Pisaura mirabilis (Clerck 1757) 

Austria 

Huber et al. 

1993 

Mygalomorphae 

Antrodiaetidae 

Antrodiaetus unicolor (Hentz 1841) 

Virginia 

AY241258 


Cyrtaucheniidae 

Aptostichus atomarius Simon 1891 

California 

AY241254 


Cyrtaucheniidae 

Aptostichus simus Chamberlin 1917 (3 
populations sampled, Bond et al. 
2001) 

California 

AF307969, 

AF307964, 

AF307960 

Cyrtaucheniidae 

Aptostichus sp. 

California 

AY241255 


Cyrtaucheniidae 

Entychides arizonicus Gertsch & Wal¬ 
lace 1936 

Arizona 

AY241257 


Cyrtaucheniidae 

Promyrmekiaphila gertschi Schenkel 
1950 

California 

AY241256 



Gutell et al. 1993) were based on extensive 
surveys of sequences as well as studies of po¬ 
sitional covariance; these are thus considered 
to be more accurate and have been used in 
recent studies modeling secondary structure in 
arthropods (e.g. Buckley et al. 2000). Existing 
secondary structure models of 16S in arach¬ 
nids are limited to a tick (Black & Piesman 
1994), which used the Clary & Wolstenholme 
(1985) model and two spider species of the 
infraorder Araneomorphae (Huber et al. 1993; 
Masta 2000), which used the Gutell model 
(Gutell & Fox 1988). To obtain a more com¬ 
plete picture of arachnid 16S secondary struc¬ 
ture, this study will examine the secondary 
structure of the other “primitive” infraorder 
of spiders, the Mygalomorphae (Coddington 
& Levi 1991), and, to a lesser extent, primi¬ 
tive liphistiid spiders, Acari (ticks) and Scor- 
piones using the Gutell model. Spider taxa 
from disparate groups were examined as part 
of a concerted effort to sample across all of 
the major clades. Table 1 summarizes the taxa 


used in this analysis, the classificatory place¬ 
ment of these taxa, and the source of their 
sequences. Decisions regarding taxon choice 
within the Mygalomorphae reflect an attempt 
to examine taxa across a number of phyloge¬ 
netic levels within this clade. 

The major objectives of this study are to: 
1) construct a secondary structure model for 
mygalornorph spiders using the preferred Gu¬ 
tell model, 2) examine trends in secondary 
structure evolution in spiders and other arach¬ 
nids, 3) analyze the pattern of molecular var¬ 
iation with respect to structure in mygalo- 
morphs and araneomorphs, and 4) assess the 
effects that differential weighting of molecular 
characters based on secondary structure has 
on phylogeny reconstruction in spiders. This 
is the first study of arachnids to examine the 
rates of variation in 16S rRNA relative to sec¬ 
ondary structure in a rigorous statistical man¬ 
ner and to analyze secondary structure trends 
across the class. It is also the first to consider 
the implications of secondary structure in se- 
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quence alignment and phylogenetic recon¬ 
struction in arachnids. 

METHODS 

DNA extractions. —Total genomic DNA 
was extracted from leg tissue using a Pure- 
gene DNA extraction kit, which comprises a 
lysis step in which ground tissue is incubated 
in Tris-EDTA buffer with SDS and Proteinase 
K for three hours, a protein precipitation step 
using potassium acetate, followed by DNA 
precipitation in isopropanol, and a 70% etha¬ 
nol wash. DNA was resuspended in Tris- 
EDTA buffer and diluted 1:100 for subsequent 
reactions. 

Mitochondrial gene PCR and sequenc¬ 
ing. —The polymerase chain reaction (PCR) 
was used to amplify the 16S rRNA gene with 
the 16S universal primers 16sar-5' (5'- 
CGCCTGTTTATCAAAAACAT-3') and 
16sbr-3' (5'-CCGGTCTGAACTCAGAT- 
CACGT-3') (Hillis et al. 1996). The primers 
16sar-5' and 16sbr-3' correspond to Drosoph¬ 
ila melanogaster mitochondrial genome posi¬ 
tions 13398 and 12887 respectively. Standard 
PCR reactions were carried out in 50 jxl vol¬ 
umes and run for 35 cycles, each consisting 
of a 30 sec denaturation at 95 °C, 30 sec an¬ 
nealing at 50°C and 45 sec ( + 3 sec/cycle) ex¬ 
tension at 72°C, with an initial denaturation 
step of 95°C for 2.5 min and a final extension 
step of 72°C for 10 min. Amplification prod¬ 
ucts were electrophoresed on a 0.8% agarose 
gel, excised from the gel and purified using 
Qiagen QIAquick gel extraction columns. Pu¬ 
rified products were sequenced from both di¬ 
rections with an ABI PRISM 377 and 310 au¬ 
tomated sequencers using the ABI PRISM 
Dye Terminator Cycle Sequencing Ready Re¬ 
action Kit with AmpliTaq DNA Polymerase, 
FS. 

Secondary Structure. —Mygalomorph 16S 
sequences generated for six taxa were crudely 
aligned (alignment of stem regions was im¬ 
proved by eye once the secondary structure 
was obtained) with sequences from eight ar- 
aneomorph taxa (Table 1) and Drosophila me¬ 
lanogaster using ClustalX (Higgins et al. 
1996). The 16S rRNA secondary structure 
was predicted for the mygalomorph Aptosti- 
chus atomarius Simon 1891 (Fig. 1) by com¬ 
parison with sequences and models from Dro¬ 
sophila yakuba (Gutell et al. 1993) and an 
araneomorph (Huber et al. 1993). We assumed 


that the location of stems and loops would be 
very similar because the LSU rRNA second¬ 
ary structure has been found to be widely con¬ 
served (Buckley et al. 2000; Masta 2000). The 
A. atomarius sequence was chosen because it 
had the highest similarity to Drosophila , fa¬ 
cilitating comparison. 

To examine evolution of secondary struc¬ 
ture in arachnids, models were also created for 
a hypervariable portion of the molecule (Fig. 
1) in distantly related arachnid taxa: the tick 
Ornithodoros moubata (Murray 1877) (Black 
& Piesman 1994), the scorpion Vaejovis car- 
olianus (Beauvois 1805) and the mesothele 
Heptathela nishihirai Haupt 1979 (Huber et 
al. 1993) (Fig. 2). All structures were drawn 
using Canvas® graphics software (Deneba 
Systems Inc.). 

Analysis of variability with respect to 
secondary structure. —Each nucleotide po¬ 
sition was assigned a single letter designating 
its structural function (S = stem, I = internal 
loop, L = external loop, C = connecting re¬ 
gion) and was coded as variable (1) or in¬ 
variable (0). In this study, stems refer to a se¬ 
ries of bonded nucleotides. Internal loops are 
those unbonded nucleotides, which occur 
within a stem; external loops occur at the end 
of stems. Connecting regions link stems. 

Variation was analyzed separately for the 
two spider infraorders, Mygalomorphae and 
Araneomorphae. For each structural category, 
variability was calculated by dividing the 
number of variable positions by the total num¬ 
ber of positions. Statistical analyses were per¬ 
formed using SAS (SAS Statistical Systems). 

To compare the pattern of variation across 
structural regions to random variation within 
the molecule, the percent variability of ran¬ 
dom blocks within the sequence was com¬ 
pared to the percent of actual variability with¬ 
in structural units (S, I, L, C). A simple 
random number generator program written for 
Mathematica (Wolfram 1996) produced nucle¬ 
otide blocks in lengths of 5 or 13 consecutive 
numbers between 1 and 503 base pairs. The 
block size corresponds to the maximum and 
minimum average sizes of the structural re¬ 
gions, and the random numbers corresponded 
to positions along the molecule. The mean 
variability in the random blocks provided the 
expected values for a Chi-Squared Test to de¬ 
termine if the pattern of variability was dis¬ 
tinct from the random model. 
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Figure 1.—Proposed secondary structure model for the 3' half of 16S of the mygalomorph species 
Aptostichus atomarius based on the Gutell et al. model (1993). Highly conserved regions indicated by 
bold-lettered nucleotides. Dashes represent Watson-Crick bonds; circles are U-G bonds. The hypervariable 
area examined in the study and modeled in Figure 3 is corresponds to nucleotide positions 237-359. 
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Figure 2.—Proposed secondary structures of a hypervariable region of 16S based on Gutell et al. model 
(1993) for 4 Arachnid taxa. Solid lines indicate tertiary interaction with strong comparative data and 
dashed lines indicate those with less support (Gutell 1996). (a) Aptostichus atomarius (Mygalomorphae) 
(b) Ryuthela nishihirai (Mesothelae) (c) tick Ornithodoros moubata (d) scorpion Vaejovis carolianus. 


Phylogenetic analysis. —Phylogenetic 
analyses were performed using PAUP* ver¬ 
sion 4.0b2a (Swofford 1999) run on a Power 
Macintosh 6500/275. The phylogenetic signal 
in the data set was evaluated using the g x sta¬ 
tistic (Hillis & Huelsenbeck 1992) based on 
100,000 random trees generated in PAUP*. 
All characters were treated as reversible, un¬ 
ordered, and all characters were initially 
weighted equally. Unambiguous gaps were 
scored as binary characters. These binary 
scorings were retained and the individual nu¬ 
cleotide positions from which they were 
scored were excluded from the analysis. Oth¬ 


erwise, gaps were treated as missing charac¬ 
ters. 

Heuristic searches were performed using 
random addition stepwise (1000 replicates) of 
taxa followed by TBR (tree bisection-recon¬ 
nection) branch swapping. Branches with a 
maximum length of zero were collapsed. Mea¬ 
sure of branch support is based on decay (Bre¬ 
mer 1988; Donoghue et al. 1992) and boot¬ 
strap analyses (Felsenstein 1985). Decay 
indices were computed using the computer 
program Autodecay (Eriksson & Wikstrom 
1996). Bootstrap values are based on 500 rep¬ 
licates using strict parsimony in PAUP*. 
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RESULTS 

Secondary structure model. —A 503 base 
pair segment of the 3' half of 16S was ob¬ 
tained for 6 mygalomorph taxa. Figure 1 
shows the proposed secondary structure for 
the mygalomorph Aptostichus atomarins 
based on the Gutell model. While the myga¬ 
lomorph model shares substantial physical 
similarity with the araneornorph and Drosoph¬ 
ila models, the underlying RNA sequences 
vary greatly except in some highly conserved 
regions (Fig. 1). The location of these con¬ 
served regions serve as anchor points for es¬ 
tablishing the molecule’s overall secondary 
structure. The structural similarity among the 
models has been maintained through compen¬ 
satory changes and by slight shifting of struc¬ 
tural regions to accommodate nucleotide sub¬ 
stitutions. 

Structural evolution in the hypervariable 
region. —Alignment of 16S sequences from 
the various arachnid taxa sampled revealed a 
highly variable length of DNA; models of this 
“hypervariable” region showed marked vari¬ 
ation in the size of loops and the length of 
stems (Fig. 2). Interestingly, bases identified 
as highly conserved (Buckley et al. 2000) are 
found adjacent to and even within this region. 
Evolution in the hypervariable segment of 
16S (Fig. 2) reveals an overall trend towards 
reduction in more advanced arachnids. We 
will describe this trend moving clockwise 
through the region (Fig. 2). The models dem¬ 
onstrate extreme variability in loop 1, which 
varies in length from 9 nucleotides in Aptos¬ 
tichus to 30 in Ornithodoros (tick). The stem 
supporting this loop is variable as well, rang¬ 
ing from 3 base pairs in Aptostichus and Hep - 
tathela to 11 base pairs in the tick. The con¬ 
necting region 2 is short in spiders as it is in 
humans and Drosophila but is longer in Vae- 
jovis (scorpion) and the tick than in spiders. 
Stem 3, which contains the sequence of four 
G’s common to all organisms, is slightly 
shorter in the spiders than in the tick or the 
scorpion. Area 4 consists of approximately the 
same number of nucleotides in the arachnids 
sampled; the morphological differences result 
from different amounts of base pairing. Re¬ 
gion 5 is maintained by two sets of bonds, 
which uphold the structural integrity of the 
stem but permit sufficient freedom for tertiary 
interactions. The stem and loop of region 6 


Table 2.-—Summary of genetic variability with¬ 
in each structural category. 


Structure 

Mygalomorphs 

Araneomorphs 

Number 

of 

posi¬ 

tions 

Per¬ 

cent 

vari¬ 

able 

Number 

of 

posi¬ 

tions 

Per¬ 

cent 

vari¬ 

able 

Stems 

65 

32% 

54 

24% 

Inner loops 

62 

45% 

62 

66% 

Loops 

103 

51% 

99 

46% 

Connecting 

217 

38% 

235 

32% 


are most reduced in Aptostichus although the 
“AUU” sequence in the loop is conserved in 
all arachnids sampled. The importance of this 
sequence is supported by Buckley et al. 
(2000), who found it to be conserved through¬ 
out insect taxa. Finally, the length of connec¬ 
tor 7 is fairly conserved, ranging from 14 nu¬ 
cleotides in Heptathela to 20 in the tick. 

Pattern of variation with respect to struc¬ 
ture. —Table 2 summarizes the pattern of var¬ 
iation within the 16S molecule for mygalo- 
morphs and araneomorphs. A Wilcoxon Rank 
Sum Test shows that the variability of each 
structural class (S, I, L, C) is statistically dif¬ 
ferent for both my gal o morphs (P < 0.05) and 
araneomorphs (P < 0.0001). 

The mean observed variability of the four 
structural classes differed from both the small 
(expected (e) = 0.29, X 2 = 40.34, P < 0.005) 
and large (e = 0.33, X 2 = 26.93, P < 0.005) 
random blocks for the araneomorphs. How¬ 
ever, the mean variability values for my gal o- 
morphs were not statistically distinct from the 
random model (small blocks: e = 0.42, X 2 = 
4.49, P > 0.10; large blocks: e = 0.43, X 2 = 
4,68, P > 0.10). 

Phylogeny reconstruction. —We consider 
the phylogenetic signal in this data set to be 
significant (gj = —0.55, P < 0.01). A strict 
parsimony analysis with all positions weight¬ 
ed equally resulted in two equally parsimo¬ 
nious trees (Figs. 3a & b), 719 steps in length 
(Cl = 0.61, RI = 0.66). Subsequent analyses 
with loop and inner loop positions first down¬ 
weighted to 0.20 and then excluded each re¬ 
sulted in one most parsimonious tree with 
slightly improved Cl and RI values (Fig. 4). 
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Figures 3 A, B.—Two equally parsimonous trees based on 16S sequences in which all positions are 
weighted equally, both 719 steps in length (Cl — 0.61, RI = 0.66). “C.” stands for the genus Cupiennius 
and “A.” for the genus Aptostichus . Locations for A. simus populations are Leo Carillo State Beach, Los 
Angeles County, CA (len). Sycamore Cove Beach, Ventura County CA (sc) and Zuma Beach County 
Park, Los Angeles County, CA (zb). 


DISCUSSION 

Secondary structure,—With the exception 
of a few hypervariable regions of the mole¬ 
cule, the mygalomorph secondary structure 
bears close resemblance to those structures 
proposed for prokaryotes (Gutell 1996; Larsen 
1992), insects (Buckley et ak 2000), and ar- 
aneomorphs (Huber et ak 1993). This supports 
the findings of Wheeler & Honeycutt (1988), 
which demonstrated that Darwinian selection 
must be operating strongly on these genes to 
maintain the functional aspects of secondary 
structure. 

The alignment of DNA sequence is com¬ 
putationally difficult (Swofford et ak 1996: 
Siowinski 1998). Comparison of the mygalo¬ 
morph model to other arthropod secondary 
structures facilitated the identification of con¬ 
served areas, which were integral to improv¬ 
ing the alignment (Buckley et al 2000). Ri- 
bosomal RNA sequences add another level of 
complexity to the problem of sequence align¬ 
ment because nucleotide site position align¬ 
ment is determined by secondary structure po¬ 
sition homology. Sequence alignments that 
account for true molecular structural homol¬ 


ogy by identifying conserved regions are most 
likely to reflect true homology and thus cor¬ 
rect gene phylogeny, a consideration that may 
not be reflected in alignment approaches that 
invoke overall similarity (e.g., Higgins et ak 
1996) or concurrent optimizations of align¬ 
ment and phylogenetic inference (e.g., Wheel¬ 
er & Gladstein 1994). 

Trends in secondary structure evolution 
in the hypervariable region,—Examination 
of the hypervariable region (Fig. 2) suggests 
a general trend toward the reduction of stems 
and loops in Araneae. The Aptostichus and 
Heptathela models exhibit a large deletion in 
the stem and loop of region 1 (Fig. 2), similar 
to that observed by Masta (2000) in Habron - 
attus oregonensis (Peckham & Peckham 
1888). The size of region 1 in the scorpion 
and the tick (39 and 53 nucleotides, respec¬ 
tively) is an intermediate between Drosophila 
(55 nucleotides) and Aptostichus (15 nucleo¬ 
tides). This suggests a pattern of evolution to¬ 
wards reduction of region 1, a possible syna- 
pomorphy of spiders. However, based on 
Wheeler & HayashPs (1998) chelicerate phy¬ 
logeny, additional 16S sequence data are re- 
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quired for at least palpigrades, amblypygids, 
and schizomids before this can be stated un¬ 
equivocally. The main exceptions to this re¬ 
ductionary trend are those areas involved in 
tertiary interactions, such as region 5 and con¬ 
necting region 7 (Fig. 2). The size of these 
areas varies little from Drosophila to humans 
to arachnids. 

The extreme variability of region 1 in spi¬ 
ders may be due to its position in the tertiary 
structure of the fully-formed ribosome. When 
addressing why a certain segment of a ribo- 
somal subunit is more variable than others, 
Simon et al. (1994) proposed that it is impor¬ 
tant to consider the position of the segment in 
the molecule. Peripheral secondary structure 
elements tend to be more variable than those 
on the interior and more frequently undergo 
modification or elimination. The hypervaria- 
bililty of this portion of 16S (Fig. 2) suggests 
that it is probably located on the periphery of 
the molecule. 

Pattern of variation with respect to struc¬ 
ture.— Within the Mygalomorphae and Ara- 
neomorphae, both the stem (S) and connecting 
(C) regions exhibited lower variability than 
the loops (I, L), although this difference was 
most distinct in the araneomorphs (Table 2). 
Because stems must base pair to maintain 
their integrity, changes that disrupt base pair¬ 
ing along the stem will be less likely to occur. 
Connecting regions are probably more con¬ 
served because they have a vital function in 
maintaining the spatial arrangement of the 
molecule. Loops are less constrained as 
changes in their sequence will not interrupt 
base pairing and at most will threaten tertiary 
interactions. Additionally, changes in length 
of the loops will have a small effect on the 
overall form of the molecule. This disparity in 
variability between stems and loops is differ¬ 
ent from what has been observed for other 
ribosomal RNA molecules. Wheeler & Ho¬ 
neycutt (1988) found stems to be more vari¬ 
able than loops for 5.8S rRNA, the small nu¬ 
clear subunit (SSU). They suggest that this 
was due to the compensatory nature of nucle¬ 
otide base changes in the stems. Mutations in 
one base pair would result in a selection pres¬ 
sure for changes in the other, causing stem 
nucleotide pairs to evolve in concert. 

Effects of differential weighting of struc¬ 
tural regions. —The results of this study dem¬ 
onstrate the importance of considering 16S 


rRNA secondary structure for phylogenetic 
reconstruction. Differentially weighting seg¬ 
ments of the 16S rRNA based on their struc¬ 
tural properties prevents the overemphasis of 
homoplasic nucleotide positions. This was 
most evident for the mygalomorph taxa. 
Down-weighting the variable loops and inner 
loops improved the phylogenetic tree, which 
recovered Aptostichus as a monophyletic 
group (47), a grouping supported by com¬ 
bined analyses of morphological and molec¬ 
ular data (Bond 1999; Bond & Opell 2002). 
This ability to fine tune a data set, omitting or 
differentially weighting hypervariable regions, 
may allow 16S sequence data to be applied 
more broadly. At the population level, highly 
variable loops would provide phylogenetic 
signal while those same loops could be dis¬ 
regarded for comparison of intermediate or 
very disparate taxa in favor of more con¬ 
served stems or connecting regions. 

Because this study utilizes only the 3' half 
of the 16S, application of our differential 
weighting scheme may require adjustment 
when considering the entire molecule. Simon 
(1991) pointed out that smaller short range 
stems (terminal stems) vary substantially 
more than long range stems (supporting 
stems) and can evolve as rapidly as some 
loops. Simon (1991) suggested that, if differ¬ 
ential weighting were employed, these two 
types of stems should be distinguished and the 
rapidly evolving shorter stems should be 
down-weighted. In the segment of 16S used 
in this study, there were 16 short range stems 
and only one strand of two long range stems. 
Therefore, the mean variability for stems 
could be considered representative of short- 
range stems. Perhaps, if the system of differ¬ 
ential weighting used in this study were em¬ 
ployed on a longer sequence of LSU rRNA, 
the two stem types would have to be distin¬ 
guished. 

Summary. —-The majority of sequence 
based phylogenetic studies rely solely on 
DNA sequence alignment algorithms and 
weighting schemes based on nucleotide com¬ 
positions. Not surprisingly, this study dem¬ 
onstrates that greater understanding of the 
processes underlying nucleotide sequence 
evolution will increase the likelihood of re¬ 
covering the correct phylogenetic relation¬ 
ships. We show that at very “deep” phylo¬ 
genetic levels gross secondary structure, the 
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Heptathela 
Pisaura 
C. getazi 
C. salei 
C. coccineus 
Dolomedes 
Pardosa 
Phoneuiria 
Clubiona 
A. atomarius 
Aptostichus sp. nov. 
A. simus (lcn) 

A. simus (sc) 

A. simus (zb) 
Promyrmekiaphila 
Entychides 
A ntrodiaetidae 


Figure 4.—Tree produced when loop and inner loop positions of 16S are down-weighted to 0.20 (520.8 
steps, Cl = 0.63, RI = 0.69) and when loops and inner loops were excluded (471 steps. Cl = 0.63, RI 
= 0.70). Bootstrap values and decay indices refer to the tree with loops downweighted. 


morphology of the molecule, may very well 
contain phylogenetic information useful in re¬ 
constructing arachnid relationships. At more 
intermediate phylogenetic levels (inter-famil¬ 
ial) we show that by down-weighting loop nu¬ 
cleotides phylogenetic signal may be im¬ 
proved. Although the conclusions drawn here 
may be limited to arachnid studies, the statis¬ 


tical approach employed may be useful for 
studies of secondary structure in other groups. 
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